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Abstract. Electron beam induced deposition with tungsten hexacarbonyl W(CO)g 
as precursors leads to granular deposits with varying compositions of tungsten, carbon 
and oxygen. Depending on the deposition conditions, the deposits are insulating 
or metallic. We employ an evolutionary algorithm to predict the crystal structures 
starting from a series of chemical compositions that were determined experimentally. 
We show that this method leads to better structures than structural relaxation based 
on guessed initial structures. We approximate the expected amorphous structures by 
reasonably large unit cells that can accommodate local structural environments that 
resemble the true amorphous structure. Our predicted structures show an insulator 
to metal transition close to the experimental composition at which this transition is 
actually observed. Our predicted structures also allow comparison to experimental 
electron diffraction patterns. 



PACS numbers: 71.15Mb, 61.46.-w, 61.43.Dq 
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1. Introduction 

Nanotechnological applications require fabrication of nanometer-sized structures on 
various substrates. Electron beam induced deposition (EBID) has emerged as a 
promising technique to make nanostructures in a size, shape and position-controlled 
manner without the use of any expensive masks [U EJ |3l IU E]- Deposits with the 
desired metal content and electronic properties can be obtained either directly by tuning 
the preparation conditions (varying electron beam energy or current or raster current) 
or by post-fabrication techniques (heating or further irradiation). Thus fabrication 
of materials with new physical and chemical properties at the nanoscale has been 
successfully achieved [3J EJ El E> U ED] • 

Transition metal carbides possess unique physical and chemical properties that 
have made them promising materials in several industrial and electronic applications. 
The composition of tungsten granular deposits obtained by decomposing W(CO) 6 as 
a precursor in the EBID process indicates that the tungsten atoms are embedded in a 
carbon (and oxygen) matrix j3]. Although investigations on the microstructure and the 
electrical transport properties have shed some light on the behavior of these systems, a 
deep microscopic understanding is still missing. 

Several theoretical studies are available on the structural and electronic properties 
of 4c? and 5d transition metal carbides [Til [12]. Nevertheless, studies on metal 
oxycarbides are scarce due to the lack of knowledge on their structures. The high 
level of carbon and oxygen concentrations up to an average of 30-40% in the EBID- 
fabricated samples indicates that a good description of the electronic structure of these 
metal oxycarbides may be obtained by suitably guessing approximate structures from 
the well-known crystal structures of tungsten carbides and tungsten oxides. This 
methodology indeed has been successful in predicting the structure of Pt 2 Si 3 derived 
from Pt2Sn3 [13]. A similar procedure for tungsten oxycarbides has been adapted 
by Suetin et ai, who investigated the structure, electronic and magnetic properties 
of some tungsten oxycarbides by constructing approximate crystal structures obtained 
from systematically replacing the carbon by oxygen in the hexagonal structure of WC 
and oxygen by carbon in the cubic structure WO3 [H] . However, a powerful evolutionary 
algorithm was recently proposed which in principle can predict the crystal structure of 
materials with any atomic composition and is not biased by the choice of initially known 
crystal structure settings JH3 ESI EI] • 

In the present work we use this evolutionary algorithm to predict structures of 
approximants to amorphous tungsten oxycarbides as obtained by the EBID process. 
By analyzing the electronic properties of our predicted structures, we find an insulator 
to metal transition at a composition close to the composition where experimentally 
such a transition has been observed [6]. We further show that the calculated X-ray 
diffraction patterns (XRD) for our structures correlate very well with the patterns 
measured experimentally [31 E] • 
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2. Method 

We approximated the amorphous tungsten oxycarbides structures obtained in the EBID 
process by large unit cells that can account for the local structural environment present 
in the experimental compositions. In order to predict these structures, we employed 
evolutionary algorithms developed by A. Oganov et al. featuring local optimization, real- 
space representation and flexible physically motivated variation operators [T5l [T6| [T7] . 
Each generation contained between 20 and 40 structures and the first generation was 
always produced randomly. Three different sets of calculation have been performed for 
each composition with differing number of initial populations and varying slightly the 
parameter (fracPerm) that controls the percentage of structures obtained by heredity 
and permutation. With all these different sets of calculation, ca. 2000 structures 
were screened for each composition. All structures were locally optimized during 
structure search using densiy functional theory (DFT) with the projector augmented 
wave (PAW) [HI HS| as implemented in VASP [121 US ED |22] ■ The generalized gradient 
approximation (GGA) in the parametrization of Perdew, Burke and Ernzerhof [23] 
was used as approximation for the exchange and correlation functional. The reported 
structures are the ones with the lowest enthalpy; the evolutionary algorithm was 
considered converged when the lowest enthalpy structure couldn't be improved during 
eight generations. We analyzed the electronic structure of the resulting structures using 
the full potential local orbital (FPLO) basis [23]. The XRD patterns were simulated by 
the Reflex module implemented in Materials Studio package. 

3. Results and Discussion 



Sample 


W 


C 


O 


approximant 


approximant 
composition 


1 


19.0 


67.1 


13.8 


w 3 c 10 o 2 


WC3.33O0.67 


2 


22.6 


56.0 


21.4 


w 2 c 5 o 2 


wc 2 . 5 o 


3 


27.5 


50.4 


22.1 


w 4 c 7 o 3 


Wd.75O0.75 


4 


31.8 


44.4 


23.8 


w 5 c 7 o 4 


Wd.400.8 


5 


34.0 


44.3 


21.7 


w 3 c 4 o 2 


WC1.33O0.67 


6 


36.9 


35.6 


27.5 


w 7 c 7 o 5 


WCO0.71 



Table 1. EBID obtained samples as reported in Ref. [5] and the corresponding 
approximant used for the structure prediction are listed. The concentrations are given 
in atomic %. The composition of the approximants normalized to the tungsten content 
is also listed. 

We first tested the method of evolutionary algorithm-based structure prediction 
using some known tungsten structures. As an example, we verified that USPEX indeed 
predicts the known hexagonal structure of WC [26]. Next, we address the problem of 
predicting crystalline tungsten oxycarbides. This has recently been discussed by Suetin 
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WC . 5 O 05 (model 2) WC . 5 O . 5 (USPEX) 




Figure 1. Predicted structures for the tungsten oxycarbide WC0.5O0.5. Left: Best 
structure obtained by relaxing guessed candidate structures. Right: USPEX result. 

et al. jH] for the examples ^NC\- x O x and ^NC^- x O x . The authors successively replace 
carbon atoms in WC with oxygen and replace oxygen atoms in WO3 by carbon atoms 
and relax the resulting structure candidates using the full potential linearized augmented 
plane wave (FPLAPW) basis set. We verified that the structure of WC0.5O0.5 with 
alternating layers of WC and WO (Fig. [I] (left)) indeed is the optimal structure also 
when relaxing different structure candidates using VASP. We then performed an USPEX 
structure prediction with the composition W2CO. This yields as optimum the structure 
shown in Fig. [I] (right). It is triclinic (PI symmetry) and it is 1.35 eV per W2CO 
unit lower in energy than the high symmetry (P — 6m2) structure obtained by relaxing 
structure candidates (Fig. [I] (left)). This indicates that indeed it is preferable to avoid 
bias by using a structure search based on evolutionary principles. 

Figs. [2] and [3] show the optimal structures we obtained. The samples with high 
carbon content show inclusions of regions that resemble diamond-like carbon (sample 1) 
or graphitic carbon (samples 2 and 4). This leads to a lower density as can be seen from 
Fig. [4] in which the density of amorphous tungsten oxycarbide approximants is plotted 
agains the metal content. On the one hand, there is a weak overall proportionality of 
density with tungsten content, illustrated by the line fitted to the six data points. On the 
other hand, samples 1, 2 and 4 fall into a lower density group (WC3.33O0.67, WC 2 5 and 
WC1.4O0.8, P from 7.8 g/cm 3 to 9.1 g/cm 3 ) which shows some phase separation between 
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WC 3.33°o.67 (sample 1) 




Figure 2. Predicted tungsten oxycarbide structures for the compositions WC3.33O0.67j 
WC 2 . 5 and WC1.75O0.75- 



low carbon content tungsten oxycarbide and regions of pure carbon, and a higher density 
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WC 1 4 O 8 (sample 4) 




wc i.33°o.67 (sample 5) 




Figure 3. Predicted tungsten oxycarbide structures for the compositions WC1.4O0.8, 
Wd.33O0.e7 and WCO0.71. 

group, samples 3, 5 and 6 (WC1.75O0.75, WC1.33O0.67 and WCO0.71, p from 10.0 g/cm 3 
to 10.7 g/cm 3 ) which is more homogeneous and more highly coordinated. 

We now investigate the electronic structure of the predicted amorphous tungsten 
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Figure 4. The density of the predicted tungsten oxycarbide structures roughly 
increases with tungsten content. 



oxycarbide approximants. We employ the FPLO basis [24J within GGA. Due to the 
large mass of tungsten, we compared scalar relativistic and fully relativistic electronic 
structure calculations. We find significant splittings due to spin-orbit coupling in all 
bandstructures, and therefore in the following we base our analysis on fully relativistic 
calculations. Fig.[5]shows a comparison of the densities of states of the six materials with 
tungsten, carbon and oxygen contributions shown in different colors. We immediately 
observe an insulator-to-metal transition between sample 2 (WC2.5O) and sample 3 
(WC1.75O0.75) which corresponds to a transition between 22% and 29% metal content. 
This is in excellent agreement with the observation of Ref. [6] where the conductivity 
measurements on the six samples (see Table [T]) showed a change from insulating to 
metallic behaviour between sample 3 and sample 4. In fact, Fig. 2 of Ref. [6] shows that 
sample 3 takes an intermediate position between clearly finite conductivity in the T — > 
limit for sample 4 and clearly vanishing conductivity in the T — > limit for sample 2. 

For comparison we present in Fig. [6] the densities of states for the predicted 
crystalline structure of our test system W2CO as well as for the known structures of 
WO2 (Ref. [25]) and WC (Ref. [26]). We observe a qualitatively different behavior for 
the three structures. 

In Figs. [7] and |8j we show the calculated bandstructures for the first four predicted 
structures. Here, we can also clearly see the transition from insulating to metallic 
behaviour upon increase of tungsten content as well as the splitting of the bands due to 
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energy (eV) 

Figure 5. Electronic density of states of the predicted tungsten oxycarbide structures. 



the spin-orbit coupling. We also observe highly dispersive bands which is a signature of 
the three-dimensionality of the systems. 

In order to check how good our simulated structures describe the amorphous 
deposits, we present in Figs. [9] and 10 the pair correlation functions of the first four 
predicted structures. The contributions of bonds involving tungsten are highlighted. 
The pair correlation function first shows carbon-carbon bonds at 1.4-1.5 A indicating 
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Figure 6. Electronic density of states of (a) the predicted structure of W2CO and of 
the known crystalline structures of (b) WO2 and (c) WC. 



that the matrix is composed of carbon atoms with sp 2 and sp 3 hybridization which is 
in accordance with the experimental evidence based on micro-Raman measurements on 
tungsten-based composites obtained in the EBID process [27J . Further analysis indicates 
that at 2.0 A, there are W-0 bonds, followed by W-C bonds at slightly larger distances. 
The first W-W bonds are seen at 2.6 A. 

Finally, we can compare our predicted structures to experiment by calculating the 



electron diffraction patterns. Fig. 11 shows the predicted electron diffraction patterns 
for electrons with an energy of 300 keV. The experimental data shown in the figure 
are from Ref. [27] . We find a very good agreement between the main peak observed 
experimentally at 4 nm" 1 and the peaks in the predicted diffraction patterns. This 
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Figure 7. Band structures of the two insulating compounds WC3.33O0.67 (top) and 
WC 2 . 5 0. 
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Figure 8. Band structures of two metallic compounds WC1.75O0.75 (top) and 
WC1.4O0.8- Compounds with higher relative tungsten content are also metallic. 
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peak corresponds to bond lengths of 2.5 A and should be related to tungsten bonds 
as the light elements contribute only insignificantly to the electron diffraction intensity 
at 300 keV. Thus, we can relate the electron diffraction pattern to the pair correlation 
functions of Figs. [9] and [10] and conclude that the W-W bonds most likely cause the 
electron diffraction peak observed experimentally. 



4. Conclusions 



By employing evolutionary algorithms we have predicted structures of approximants to 
EBID-based amorphous tungsten oxycarbides. By analyzing the electronic structure, 
pair correlation functions and diffraction patterns of our predicted structures for different 
compositions of W, C and O, we find very good agreement with the experimental 
observations; an insulator-to-metal transition is observed at a concentration close to 
the experimental concentration at which this transition is reported and we explain the 
prominent peak observed in micro- Raman spectroscopy as caused by W-W bonds. 

The use of genetic algorithms to predict and simulate amorphous nanodeposits 
opens the possibility to understand the microscopic origin of the behavior of these 
systems. This was up to now very limited and mostly restricted to phenomenological 
models. With this tool at hand, we believe that important progress can be made in this 
field. 
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Figure 9. Pair correlation functions of the compounds WC3.33O0.67 (top) and 
WC2.5O. 
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Figure 10. Pair correlation functions of the compounds WC1.75O0.75 (top) and 
Wd.4O0.s- 
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Figure 11. Calculated electron diffraction intensities (lines) for an electron energy of 
E = 300 keV, compared to the measured spectrum of Ref. [57] . 



